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The Lifshitz theory of dispersion forces leads to thermodynamic and experimental incon- 
sistencies when the role of drifting charge carriers is included in the model of the dielectric 
response. Recently modified reflection coefficients were suggested that take into account 
screening effects and diffusion currents. We demonstrate that this theoretical approach 
leads to a violation of the third law of thermodynamics (Nernst's heat theorem) for a 
wide class of materials and is excluded by the data from two recent experiments. The 
physical reason for its failure is explained by the violation of thermal equilibrium, which 
is the fundamental applicability condition of the Lifshitz theory, in the presence of drift 
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and diffusion currents. 
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1. Introduction 

In the last few years the Casimir effect has attracted increasing attention due to 

prospective apphcations in both fundamental physics and nanotechnology. Sixty 

years ago H. B. G, CasimiJUmade his famous discovery that two electrically neutral 

parallel ideal metal plates spaced at some separation a in vacuum attract each 

other. Casimir explained this effect as due to the alterations in the spectrum of 

zero-point oscillations of the electromagnetic field introduced by the presence of the 

plates. In the simplest approach, the Casimir effect can be described theoretically 

using quantum field theory with boundary conditions. The case of real material 

2 

plates was considered by Lifshitz, who described material properties by means of 
a frequency-dependent dielectric permittivity. The Lifshitz theory was successfully 
applied to the interpretation of measurement data in several experiments measuring 
the Casimir force.^^ 

Further investigations revealed that the Lifshitz theory at nonzero temperature 
leads to problems when the relaxation of free charge carriers is included in the 
model of the dielectric response. This was shown to lead to both thermodynam- 
ically and experimental inconsistencies.^"^ Problems arise in the zero-frequency 
contribution of the Lifshitz formula when the free charge carriers are described by 
the dielectric permittivity of the Drude model. The Lifshitz theory was found to be 
thermodynamically consistent and in agreement with the experimental data if the 
free el ectrons in metals are described by the dielectric permittivity of the plasma 
j-i^Q(jgj[T0]-[12j Pqj, (dielectric and semiconductor materials, consistency with thermo- 
dynamics and experiment is achieved if the dc conductivity is neglected or charge 
carriers are described by means of the plasma model depending on whether the 
concentration of charge carriers is below or above the critical value, respectively.^!^ 

Recently, an alternative a ppr oach to the description of free charge carriers in the 
Lifshitz theory was suggestedl^'"'!^ which takes into account screening effects and 
diffusion currents. Within this approach, the macroscopic characteristic of the plate 
material (the dielectric permittivity) is supplemented with a microscopic quantity 
(the density of free charge carriers) . Below we consider the most typical configura- 
tion of two thick parallel plates (semispaces) and demonstrate that this approach 
is thermodynamically and experimentally inconsistent. The reason for this failure 
is the violation of thermal equilibrium which is the basic applicability condition of 
the Lifshitz theory. 

The structure of the paper is as follows. In Sec. 2 we briefly formulate problems 
which arise when the Lifshitz theory is applied to materials with nonzero conduc- 
tivity. Section 3 contains the formulation of approaches taking into account the 
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charge screening of free carriers. In Sec. 4 it is demonstrated that these approaches 
are in conflict with thermodynamics. Section 5 is devoted to the comparison of 
theoretical resuhs obtained with the inclusion of the screening effects and available 
experimental data. In Sec. 6 the reader will flnd our conclusions and discussion. 

2. Problems of the Lifshitz Theory in AppHcation to Materials 
with Nonzero Conductivity 

At nonzero temperature all materials have a nonzero conductivity. For metals and 
metallic-type semiconductors this conductivity can be rather large, and does not 
go to zero when the temperature vanishes. For dielectrics and some semiconductors 
(intrinsic ones and those with dopant concentration below critical) conductivity is 
much smaller than for metals, and goes to zero together with temperature. In the 
Lifshitz theory, the free energy per unit area in the configuration of two semispaces 
described by the dielectric permittivity £{uj) is given by 

ln[l-r^E(iei,fci)e-2-9<]}. 



1=0 



Here, a is the separation distance between the semispaces, is the magnitude of 
the wave vector in the plane of boundary plates, the primed sum adds a multiple 
1/2 to the term with ^ = 0, = 27rfcsTZ/fi, with / = 0, 1, 2, . . . are the Matsubara 
frequencies, and fc^ is the Boltzmann constant. Equation ([T]) is derived under the 
condition that the plates (semispaces) are at a temperature T in thermal equilibrium 
with the environment. The reflection coefficients for two independent polarizations 
of the electromagnetic field coincide with the Fresnel ones calculated along the 
imaginary frequency axis, 

rTM(Ki,k_L) = —, rTE[iU,k±) = —, (2) 

£iqi +ki qi + ki 



where 



qf = kl + ^, kf^kl+si^, £, = e(iCO- (3) 



Originally, Eqs. H])-® were mostly applied to dielectrics with the dc conduc- 
tivity neglected. For such materials the dielectric permittivity can be written in the 
fornJl3 

where LUj ^ are the oscillator frequencies, fj are the oscillator strengths and 
are the relaxation parameters. This equation describes the dielectric response of 
core electrons. The free energy ([T]) with the dielectric permittivity ([4]) is in perfect 
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agreement with thermodynamics. Specifically, the Casimir entropy vanishes with 
temperature!!^ 

Problems with thermodynamics arise when the conductivity a of the plate mate- 
rial is taken into account. For metals and metallic- type semiconductors the dielectric 
permittivity can be modelled by means of the Drude model, 

m = e(io + = em + (5) 

where ujp is the plasma frequency and 7 is the relaxation p aram eter of free electrons. 
Simple expressions for the quantities entering Eq. ([5]) areff^ 



aiC ujI^ , a(0) = M|en, 6 

1 + i- m 

1 

where cr(0) is the dc conductivity, e and m are the charge and effective mass of the 
electron, n is the density of charge carriers and is their mobility. The substitu- 
tion of the dielectric permittivity ([5]) into Eqs. ([I])-® results in a negative Casimir 
entropy at T = K depending on a for metals with perfect crystal lattices,^^ which 
is in violation of the third law of thermodynamics (the Nernst heat theorem) . Re- 
cently it was claimed!^ that Ref. HI is wrong because it uses the theory of the 
normal skin effect for metals with perfect crystal lattices at T — > 0, while in this 
situation one must use the theory of the anomalous skin effect. This objection is, 
however, incorrect. With the decrease of T the application region of the normal skin 
effect does become narrower and the application region of the anomalous skin effect 
widens. However, at any T > 0, there exists a frequency region near zero frequency, 
where the normal skin effect is applicable. Keeping in mind that the violation of 
the Nernst heat theorem originates entirely from the zero-frequency term of the 
Lifshitz formula, one concludes that the permittivity of the normal skin effect is 
appropriate for the evaluation of this term at low T when conside ring th e thermo- 
dynamic consistency of the Lifshitz theory. There is a sugg estioi|2QEl] 

to satisfy 

the Ner nst he at theorem by the inclusion of impurities, but this does not solve the 
problen|22123] because according to quantum statistical physics for perfect crystal 
lattices the Casimir entropy at zero temperature must be equal to zero. 

Another problem is the experimental inconsistency of the Lifshitz theory com- 
bined with the dielectric permittivity ([5|) . This was demonstrated in a series of three 
successive experiments on th e dynamic determination of the Casimir pressure be- 
tween two gold plates.^ ^ l ^^ l For example, in Fig. 1 we demonstrate t he c omparison 
between the measurement data of the most precise third experimenlPEl and theo- 
retical Casimir pressures computed using the Drude model ([5]) and the generalized 
plasma- like model of Refs. [T0Ul2l [dielectric permittivity ^ with 7 = 0]. In Fig. 1(a) 
the Casimir pressures computed using the Drude model approach and the general- 
ized plasma-like model are shown as dark-gray and light-gray bands, respectively. 
The experimental data are shown as crosses. The widths of the bands and the sizes 
of the arms of the crosses are determined at a 95% confidence level. In Fig. 1(b), the 
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P (mPa) P"i"r - P'^P' (mPa) 




Fig. 1. (a) The crosses show the measured mean Casimir pressures together with the absolute 
errors in the separation and pressure as a function of the separation. The theoretical Casimir 
pressures computed using the generalized plasma-like model and the optical data extrapolated by 
the Drude model are shown by the light-gray and dark-gray bands, respectively, (b) The differences 
of the theoretical and the mean experimental Casimir pressures between two Au plates versus 
separation are shown as dots. The theoretical results are calculated using the Lifshitz theory at 
room temparature using the generalized plasma-like model (the dots labeled 1) and the Drude 
model approach (the dots labeled 2). The solid and dashed lines indicate the boundaries of 95% 
and 99.9% confidence intervals, respectively. 



differences of the theoretical Casimir pressures computed using Eq. (O, P||^®°'"(a), 
and the mean experimental Casimir pressures, P'^''P*(a), are shown as dots labeled 2. 
The differences of the theoretical pressures computed using the generalized plasma- 
like permittivity [Eq. ([5|) with 7 = 0], Pg^°°'^{a), and P^^^*-{a) are shown as dots 
labeled 1. The solid and dashed lines indicate the boundaries of 95% and 99.9% 
confidence intervals, respectively. As can be seen in Fig. l(a,b), the generalized 
plasma-like model is consistent with data, whereas the Drude permittivity taking 
into account the relaxation properties connected with a drift current of conduction 
electrons is experimentally excluded. 
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Fig. 2. (a) The measured force differences are shown as crosses versus separation a. The solid 
and dashed lines present the theoretical results calculated with the neglected and included dc 
conductivity of high-resistivity Si in the dark phase, respectively, (b) Theoretical minus mean 
experimental differences of the Casimir force for the Lifshitz theory with neglected (label 1) and 
included (label 2) dc conductivity of high-resistivity Si in the dark phase. The solid lines indicate 
95% confidence intervals. 
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Now we consider dielectric plates, whose conductivity vanishes with T, and take 
into account the small dc conductivity using the first equality in Eq. ([5]) in the 
calculation of the Casimir free energy and entropy [the second equality in Eq. ([5]) is 
related only to metals]. In this case the Casimir entropy at T = takes a positive 
value depending on a,5 13 25 21 i 

the Nernst heat theorem is violated. Theo- 
retical results computed with inclusion of the dc conductivity of the dielectric are 
in disagreement with the experimental results on the measurement of the C asimir 
force between a gold sphere and a silicon plate illuminated with laser pulses .1221281 
As an illustration, in Fig. 2(a) the mean measured differences of the Casimir forces 
in the presence and in the absence of light on the plate (the absorbed power is equal 
to 4.7mW), -Fj^fj*, are shown as crosses (the arm sizes are determined at a 95% 
confidence level). The theoretical differences, F^fg™, shown by the solid line, are 
computed by using the dielectric permittivity of the generalized plasma-like model 
in the presence of light and neglecting the dc conductivity in the absence of light 
on the plate. The theoretical differences shown by the dashed line are computed 
by using the dielectric permittivity ([5]) with the respective values of a in the pres- 
ence, and in the absence, of light, i.e., with included dc conductivity of the silicon 
plate also in the dark phase. As is seen in Fig. 2(a), the theory taking into account 
the conductivity of Si in the dark phase using Eq. ([5]) is experimentally excluded, 
whereas the theory neglecting the dc conductivity of dielectric silicon is consistent 
with the data. 

The comparison of experiment with different theories can also be done in an- 
other manner. In Fig. 2(b) we plot as dots the theoretical minus mean experimental 
differences of the Casimir forces, i^^iff"' ~ ^dis''' presence and in the absence 

of light. For the dots labeled 2 the values of i^^iff ""^ ^^'^ computed with the dc con- 
ductivity of Si included in the dark phase, and for dots labeled 1 the dc conductivity 
is neglected. The solid lines indicate the boundaries of 95% confidence intervals. As 
is seen in Fig. 2(b), the theory including the dc conductivity in the dark phase is 
excluded by the data over a wide separation region, whereas the theory neglecting 
the dc conductivity is experimentally consistent. 

The physical reasons why the use of the dielectric permittivity ([5]) in the Lif- 
shitz theory leads to problems with thermodynamics and experiment are discussed 
in Ref. [29l The point is that this permittivity describes the influence of the drift 
current of conduction electrons, an irreversible process. It is accompanied by Joule 
losses and the heating of the Casimir plates. To keep the temperature constant, 
one should allow the existence of an unidirectional flux of heat from the plates to 
the heat reservoir This is a situation out of thermal equilibrium.SH Thus, the 
dielectric permittivity ([5|) describes processes violating thermal equilibrium which 
is the basic applicability condition of the Lifshitz theory. It is thus not surprising 
that in combination with the permittivity ((5|) , the Lifshitz theory becomes thermo- 
dynamicly and experimentally inconsistent. 
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3. Attempts to Generalize the Lifshitz Theory Through the 
Inclusion of Screening Effects 

As discussed in Sec. 2, to avoid problems in tlie application of the Lifshitz theory 
to real materials one should neglect the dc conductivity of dielectrics and describe 
charge carriers in metals using the generalized plasma-like dielectric permittivity. 
Another approach to the resolution of problems arising in the Lifshitz theory when 
it is applied to materials with nonzero conductivity is by including the effect of 
Coulomb screening by mobile charges. As discussed above, at nonzero temperature 
all materials contain charge carriers which can move through the medium. In the 
presence of mobile charge carriers the standard Coulomb potential of some external 
charge is replaced with the so-c alle d screened Coulomb potential having a Yukawa- 
type dependence on separation.I32] The effect of screening leads to penetration of a 
static electric field into materials with nonzero conductivity to a depth of the screen- 
ing length which is very small for good metals (of order the interatomic distance), 
but can be much larger for semiconductors. 

The response of charge carriers to a spatially variable, sinusoidally varying, 
external electric field taking into account scattering can be described by Boltzmann 
transport equation. This approach was recently use to describe the interaction of 
the fluctuating electromagnetic field with mobile charge carriers in a semiconductor 
plate and to find the modified reflection coefficients for the TM and TE modes. 
The resulting coefficients were substituted into the standard Lifshitz formula ([1]). 
Here, we consider the applicability of this approach to different materials leaving the 
problems of its consistency with thermodynamics and experiment for the following 
sections. 

The Boltzmann transport equation is used for the description of irreversible 
nonequilibrium processes, and is not symmetric with respect to time reversal. As 
a result, it describes processes which lead to an increase of entropy.!^ In the case 
under consideration this equation includes both the drift current of charge carriers 
j and the diffusion current eDWn, where D is the diffusion coefficient. As a result, 
the modified TM reflection coefficient takes the for nP 

em + ki + kpji £, (£/ - ei) 

Here, ei and qi are deflned in Eq. ([3]), ii = e{i£_i) is defined in Eq. ki is given by 
Eq. ^ where is replaced with ii and 



m = jyte) = 



2 £o £i 

K 



£l £l — £l 



1/2 

£o = £(0). (8) 



The screening length 1/k in Eq. ^ is different for different types of statistics of 
charge carriers (see below). The modified TE reflection coefficient, r™^'^{i£^i, k±), is 
given by the standard expressions ([2]) and ([3]) , where the permittivity ei is replaced 
by 
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It has been showiJ33l that the problems of the standard Lifshitz theory discussed 
in Sec. 2 are mathematically connected with the discontinuity of the zero-frequency 
term of the Lifshitz formula with respect to frequency and temperature at the point 
(0,0). The modified reflection coefficients at zero frequency take the form 



r'^^^{0,k^) = ^-^ :±Z:L :\ r??^<^(0,fc^) = 0. (9) 



Note that the TM coefficient in ([9|) was first obtained by PitaevskiP^ within the 
theoretical approach taking into account the penetration of the static component 
of the fluctuating field in the material of a wall in the atom-wall (Casimir-Polder) 
interaction. According to the modified TE reficction coefficient at zero frequency 
takes the same value as in the Drude model approach. Mathematically, this suggests, 
that for metals the theoretical approach taking into account screening effects should 
face the same difficulties as the Drude model approach (see Sec. 2) . 

The reflection coefficients r^MTsi^T^^) '^^^ obtainec j^^ l ^^ l as the standard 
Fresnel reflection coefficients for an uniaxial crystal, 

U 1 \ ^/^Xl^zlll — kzl U f-f 1 \ Ql ~ kxl 



^j£xi£zm + kzl qi + kxl 



where 



£2 £2 

kll^k\+ezl^, Exl^ Exilil), Ezl = £z{\S,l), (11) 



if one introduces the dielectric permittivity depending on the wave vector k 



-1-7 



^2 



£x0=£x{Q) = £0, £zO = £z{0) = eo\l + -^j. (12) 

This means that the standard Lifshitz formula (at least its zero-frequency term) is 
applied to spatially dispersive materials. Such an application is controversial and has 
been debated in the literature.^^Slxhe modified refiection coefficients 7'TM^TE(i6i ^i) 
at any frequency a re a lso claimed to be obtainable in terms of spatially nonlocal 
dielectric functions.!!^ It is easily seen, however, that there are no such dielectric 
functions £xi and £zi which transform the reflection coefficients (fTO| for an uniaxial 
crystal into »^TM^TE(i^^ ^-i-)- To prove this, we first require that ^^^eIi^j ^-l) be equal 
to r™^{\^i, k±). For this purpose one must put £xi = Si- Then, keeping in mind that 
ii goes to infinity when ^ — > 0, one finds that in this limiting case r^^yi{i-Ch k±) in 
Eq. (|10p goes to unity independent of the functional form of £zi- Thus, there is no 
such function EzI that would lead to the required equality r^^ui^^ k±) — »'tm^(0, k±), 
as defined in Eq. 

In Ref . [36l the modified reflection coefficients r^^MTEi^^i^ k±) at any frequency 
are expressed in terms of two dielectric functions ^^^(i^: k±) and £^(1^;, k±) defined 
in the random phase approximation.- The obtained expression for £z(i^z,fc_L) is, 
however, incorrect. It does not transform into the dielectric permittivity of dielectric 
material £(i^i) defined in Eq. ([4]) in the limiting case n — > 0. This is caused by several 
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mistakes in sign made in Ref . f36| and by the disagreement of phase multiples in the 
TM reflection coefficients used by the authors of Ref. |36] and in the method of 
random phase approximation. After correcting these errors one obtains 



Here, fc; is defined by Eq. ([3]) where ei is replaced with £x{^^i, k±). The latter coin- 
cides with e(i^/) from Eq. (O, i.e., it does not depend on k±. It is easily seen that 
in the limiting case n — > the permittivity (|13p transforms into ei = e{i^i). 

We emphasize that the coefficient r'^^{i^i, k±) obtained in the random phase 
approximation does not coincide with the coefficient r^yi{iS,i, k±) for an uniaxial 
crystal, as given by Eq. p^ . It is notable also that the representation for the 
reflection coefficient r^yii^^i' ^J-) t>y means of the dielectric permittivity ez(i6, k±) 
at both zero and nonzero frequency is a phenomenological one. The point is that in 
the presence of a gap between the plates the translational invariance in z-direction, 
i.e., perpendicular to the plates, is violated and ez(i^;, fcj.) is an ill-defined quantity. 

Phenomenological reflection coefficients at any frequency, depending on two spa- 
tially nonlocal dielectric permittivities Sxi, £zi and the scree ning length 1/k, were 
also suggested without refer to Boltzmann transport equation.^1^ At zero frequency 
they take the same values as in Eq. ([9]). This approach is a crude approximation. 
Given the absence of translational invariance along the z-axis mentioned above, it is 
impossible to define the dielectric permittivity Szi depending on the frequency and 
the wave vector. Additionally, specular reflection of charge carriers on the boundary 
planes was assumed. However, for spatially dispersive materials the scattering of 
carriers is neither specular nor diffuse!^ Because of this, it was concludecPS] that 
this approach does not contain self-consistent checks of its accuracy, and hence could 
only be justifled based on agreement with the experimental data and fundamental 
physical principles. In the next two sections we demonstrate, however, that theoret- 
ical approaches taking the screening effects into account are in disagreement with 
thermodynamics and available experimental data. 

In the end of this section we discuss the claimed application region of the mod- 
ified reflection coefficients taking into account screening effects. This is a familiar 
subject when describing conducting materials in an external electromagnetic field 
(we recall that Ref. [IB, considers the fluctuating fleld in a similar way as an external 
one). In Ref. [151 the Boltzmann transport equation is applied in the nondegenerate 
continuum limit for sufficiently low density of charge carriers n (intrinsic semi- 
conductors). In this case the screening length is given by a specific Debye-Hiickel 
expression. 




(13) 



1 



1 




(14) 
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In fact this expression is applicable to charge carriers obeying Maxwell-Boltzmann 
statistics. It is obtained from the general representation for the screening lengthp^ 



^"^ (15) 



K Y 47rcr(0) ' 

if one uses <t(0) from Eq. ^ and Einstein's r elatioiil^l^ D / jjL ^ kgT/lel valid in 
the case of Maxwell-Boltzmann statistics. 

However, the application region of the modified reflection coefficients r™^r^^ 
with the Debye-Hiickel screening length (|14|) is not restricted to only intrinsic semi- 
conductors, but they are applicable to all materials where n is not too large so that 
charge carriers are described by Maxwell-Boltzmann statistics. This means that it 
is legitimate to apply the coefficients ?'tm^te doped semiconductors with dopant 
concentration below critical, and to solids with ionic conductivity, etc. 

Reference 15 claims that its approach "does not apply to metals, where the 
electron density is sufficiently large and the electron gas is degenerate" . The Boltz- 
mann transport equation, however, is equally applicable to classical and quantum 
systems. The only difference one should take into account is the type of statis- 
tics. In metals and metallic-type semiconductors charge carriers obey the quantum 
Fermi-Dirac st atistics . Substituting Einstein's relation valid in the case of Fermi- 
Dirac statistics jlI^ISH/)/^ = 2£'i?/(3|e|), where Ep — fiujp is the Fermi energy, into 
Eq. one arrives at the Thomas- Fermi screening lengthjS^l 



11/ £oEf 
K ktf V DTre^n 

With this definition of the parameter k, it is justified to apply the modified reflection 
coefficients rJJJ^'^rpg to metals. Thus, with the proper definition of the screening 
length, the suggested generalization of the Lifshitz theory, if at all meaningful, 
should be applicable to any material. The problem remains whether the screening 
effects are relevant to the Casimir force. This is discussed in the next sections. 



4. Inclusion of Screening Effects and Thermodynamics 

Here we perform the thermodynamic test of the Lifshitz formula ([T]) combined 
with the modified reflection coefficients from Eqs. ([7]), ([8]) and r^^'^. For this 
purpose the Casimir entropy at zero temperature will be calculated analytically. It 
is convenient to introduce the dimensionless variables y = 2aqi and Q = ^i/lOc = 
2aS,i I c. Then the modified Casimir free energy takes the form 



87ra2 



(=0 



(17) 
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In terms of dimensionless variables the reflection coefBcient ^ can be rearranged 
as 



'^TM(iO,yj = — — (lo) 



where 



1/2 

, Ka = 2aK. (19) 



£;(e; - £i) 

Note that aU dielectric permittivities here are functions of iuJcCi- Below we do not 
use the explicit expression for the reflection coefficient fTE(iC/j2/) because it coin- 
cides with the standard one, as defined in the Drude model approach, which was 
considered in detail in Ref. 141 

Let us determine the behavior of the Casimir free energy (|17p at low tempera- 
tures. We begin with the case of metals where k = ktf- For all metals the screening 
length is very small. As a result, at any reasonable separation distance between the 
plates, the dimensionless parameter Ka defined in (jl9p is very large and the inverse 
quantity /3a = '^1 i^a ^ 1 can be used as a small parameter. Expanding the reflection 
coefficient (|18p up to the first power of the parameter (3a one obtains 

^TM (iO, 2/) = rTM(iO, y) - 2/3a Zi + 0{Pl), (20) 



Zi 



lei{ei-eif y{y' - Cf) 



£o£i [eiy + VF + W^lK/ 



212 • 



In 



where f tm (iO ; y) is the standard TM reflection coefficient calculated with the di- 
electric permittivity e{iu>cCi)- [It is given by Eq. (|18p with the third terms in both 
the numerator and the denominator omitted.] From Eq. ()20|) one arrives at 

- r^^^'\iCi,y) e-y] = In [l - f'^M,y) + ^Pa J^^^^f'f/\ + 0{Pl). 

(21) 

Now we substitute (|2ip and the respective known expression for the TE 
contributiorPl into ([TT]). Calculating the sum with the help of the Abel-Plana for- 
mula, we obtain in perfect analogy to Ref. [H 

r^^°\a, T) = T,,{a, T) - M_ j^y dy In [l - ^^,,^(0, y) e'l 

+ T^'<\a,T)+PaT^f^\a,T), (22) 

where J-'^''\a^ T) is determined by Eq. (17) in Ref. 4, It goes to zero together with its 
derivative with respect to temperature, when T — > 0. The quantity JF'^'^^(a, T) origi- 
nates from the second contributions on the right-hand sides of ((20|) , ((2T|) . Using the 
Abel-Plana formula,"^ it can be easily seen that T'^'^\a, T) = E'^'^\a) + 0{T^/T^ff) 
at low T. The Casimir free energy !Fgp{a,T) is defined by substituting the dielec- 
tric permittivity of the generalized plasma-like model [Eq. ([5]) with 7 = 0] into the 
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Lifshitz formula. It was found in Refs. fTTl[T2l and the respective thermal correction 
was shown to be of order (T/Tcs)^ when T ^ 0. The TE reflection coefficient cal- 
culated using the generalized plasma-like model at zero frequency entering (j22p is 
given by 



cy - ^/la^u;; + c^y"^ 



rTE,<,p(0, y) = \ (23) 



cy -h A/4a'w; + c2y2 



Finally, calculating the Casimir entropy with inclusion of the screening effects 

"°d(, 



S (a,T) = — , (24) 



using Eq. in the limit T ^ 0, one obtains 

7 />00 

5'"°'^(«,0) = ydy ln[l-4E,gp(0,y)e-^] < 0. (25) 

Thus, the Nernst heat theorem is violated and the theoretical approach leading 
to the modified reflection coefficients is thermodynamicly inconsistent. Note that 
this result is obtained for metals with perfect crystal lattices. In the presence of 
impurities the Casimir entropy abruptly jumps to zerd^ at T < 10~^K. This, 
however, does not solve the problem of consistency with quantum statistical physics 
as discussed in Sec. 2. 

Next we consider the low-temperature behavior of the Casimir free energy for 
dielectric materials which includes screening effects. This is also relevant for semi- 
conductors with concentration of charge carriers below the critical value. For these 
materials n is relatively small and k = kdh, as defined in Eq. ([T4|) . should be used. 
The derivation of the low-temperature behavior can be performed_^as in the case 
of two dielectric semispaces with the inclusion of dc conductivity.-^ For dielectric 
materials the small parameter /3; is given by 

A = ^ (^>l), -(O)-exp(-^), (26) 

where C is some constant having a different meaning for different classes of di- 
electrics. The parameter /?; goes to zero when the temperature vanishes. Then we 
expand the modified reflection coefficients J'xm^te i^O j v) with Z > 1 in powers of the 
small parameter /3;, 

mod/-/- \ \ , a + (£( — 2)Cj^] q'2\ 

Ttm {<uy) ^ r^MiKuy) + I3i^= — —^=====-+0{l3i), 

yjy^ + {ei - i)(:t[£iy + vy + ('^i - 1)C;1^ 

y[y - \/y'^ + (e; - '^W 
^y^ + {si-i)Ct[y + ^y^ + {ei-i)(i 

Here, the reflection coefficients ^tm.te are defined by Eq. ^ with the dielectric 
permittivity of core electrons Q. The Casimir free energy J-{a,T) calculated with 
the coefficients ttm.te vanishes with temperature as ~ (see Ref. [5]). 



r^^^'iKuy) = rTE(io,y)+A .rx ' .^'^r 
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Now we substitute ([77)) in Eq. (fT7|) and arrive at the following expression for the 
Casimir free energy taking the screening effects into account: 



y dy In 



mod'^ 
'TM 



[0,y)e~y\ +Li3(r2) + Q(T)j, 

(28) 

where Li„(z) is the polylogarithm function, Q{T) vanishes exponentiallj0 when 
T 0, ro = (eo - l)/(eo + 1) and 



r-TM(0,y) 



eov^y^ + (2aKDH)^ - y 



(29) 



£Q\/y'^ + (2aKDH)^ + y 

is the reflection coefficient ^ expressed in terms of the dimensionless variables. 
Calculating the negative derivative of both sides of ([28]) with respect to T, we 
obtain the asymptotic behavior of the Casimir entropy at low temperature 



Sa^eoT 




ev - r^^\0, y) Vv^ + (SaKDH)^ [eq ^2/^ + C^ai^m)^ + y? 



(30) 



where S{a, T) is defined using J-{a, T) and, thus, vanishes when T — > 0. It is easily 
seen that the last three terms in curly brackets on the right-hand side of this equation 
go to zero when T goes to zero for any dielectric material. 

The behavior of the first two terms in the curly brackets on the right-hand side of 
pop when T goes to zero is more involved. If n{T) exponentially decays to zero with 
vanishing temperature (as is true for pure insulators and intrinsic semiconductors), 
then according to ([14]) so does kdh- As a result, r^^{Q, y) tq and the first two 
terms in the curly brackets cancel. Then the Casimir entropy 5'"°'*(a,r) goes to 
zero when T vanishes following the same law as S{a, T), i.e., as T^. This means that 
for insulators and intrinsic semiconductors the formalism under consideration is in 
agreement with the Nernst heat theorem. 

However, there is a wide class of dielectric materials (such as doped semicon- 
ductors with dopant concentration below critical, dielectric like semimetals, certain 
amorphous semiconductors, and solids with ionic conductivity) for which n does 
not go to zero when T goes to zero. Although a{0) goes to zero exponentially fast 
for all dielectrics when T goes to zero, for most of them this happens due to the 
vanishing mobility [see Eq. ([6])]. For instance, the conductivity of Si02 discussed in 
Ref. 14 is ionic in nature and is determined by the concentration of impuritues. For 
all such materials, in accordance with Eq. (|14p . kdh ^ oo when T ^ 0. As a result, 



' TM 



1^(0, y) ^ 1 when T goes to zero in accordance with (|29p . In this case we obtain 
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from ([50)1 



S""°'^(a,0) 



[C(3)-Li3(r2)] >0, 



(31) 



i.e., the Casimir entropy is positive and depends on separation distance in violation 
of the Nernst heat theorem. This means that for a wide class of dielectric materials 
the proposed approa ch t aking the screening effects into account is thermodynam- 



dc conductivity included. 

We emphasize that the existence of dielectric materials for which n does not 
go to zero but /x does go to zero when T vanishes demonstrates that the reflection 
coefficient ((T8)) at ^ = is ambiguous. In reality, for such materials r^^i^^ k±) ^ 1 
when T and simultaneously vanish. This is because kdh oo when T — > 
in disagreement with physical intuition that there should be no screening at zero 
mobility. This ambiguity is connected with the break of continuity of the reflection 
coefficient f™^'^(i^, fcj^) at the point ^ = 0, T = 0. If one takes the limit T ^ first, 
keeping ^ = const ^ 0, the standard Fresnel reflection coefficients ttm from Eq. ([2|) 
with no screening are reproduced. This property is preserved in the subsequent 
limiting transition ^ ^ 0. Thus, as was already noted above, the violation of the 
Nernst heat theorem is caused by the break of continuity of the reflection coefficients 



Recently it was claimed^ that the nonlocal approach leading to the reflection 
coefficient ([9]) with k = kdh satisfies the Nernst theorem, specifically, for solids 
with ionic conductivity which is the conductivity of activation type. To prove this, 
Ref. il6l arbitrarily separates the thermal dependence of (7(0) in Eqs ([26|) from 
the mobility fi and attributes it to the "effective density of charges, which are able 
to move" . This transfer of the temperature dependence from /x to n is incorrect!^ 
because the commonly used density of charge carriers n producing the effect of 
screening in ionic conductors is an independently measured quantity, which does 
not vanish with T. Independent measurements of all three quantities, conductivity, 
charge carrier concentration and mobility, demonstrate that "mobi lity has the dom- 
inating influence upon the conductivity-temperature dependence" Nevertheless, 
Ref . 133] expresses doubts concerning the existence of dielectric materials whose den- 
sity of charge carriers does not go to zero when T vanishes, while the conductivity 
goes to zero due to vanishing mobility. It is common knowledge, however, that for 
semimetals of the dielectric type the Fermi energy is at a band where the density of 
states is not equal to zero. The number of charge carriers (electrons) near the Fermi 
surface in such dielectric materials is fixed because it is determined by the struc- 
ture of the crystal lattice. Thus, the density of charge carriers is nonzero at any T 
including T = 0^ The same holds for certain of the amorphous semiconductors^ 
and for doped semiconductors with dopant concentration below critical. 

Recent Ref. [36] claims that it explicitly shows the satisfaction of the Nernst theo- 
rem in the Lifshitz theory with the modified reflection coefficients ?'xm*te introduced 




same way as the standard Lifshitz theory with the 
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in Ref. [151 According to Ref. [36l the Nernst theorem is satisfied "in systems with 
low density of carriers (intrinsic semiconductors, dielectrics, etc)" . However, the 
presented proof uses an assumption that "the carrier density vanishes as T — s- 0" I^SI 
Thus, wide classes of dielectric materials for which the density of charge carri- 
ers does not vanish when the temperature vanishes (doped semiconductors with 
dopant concentration below critical, semimetals of dielectric type, some amorphous 
semiconductors and solids with ionic conductivity) are simply excluded from con- 
sideration. As a result, the proof of Ref. i36iis in fact applicable to only insulators 
and intrinsic semiconductors, where the Fermi energy at T = lies in a band gap. 
For these materials the density of states at the Fermi energy is equal to zero. Thus, 
for insulators and intrinsic semiconductors the Nernst theorem is satisfied, as was 
proved above on the basis of Eq. ([30]) . However, for the majority of dielectric mate- 
rials n does not go to zero with vanishing T leading to the violation of the Nernst 
theorem in the modification of the Lifshitz theory proposed in Ref. [15] and also in 
Refs. [H M 

There is a remark in the literatur that the violation of Nernst 's theorem for 
dielectric materials where n does not go to zero with vanishing T is "a pure mis- 
understanding" because "The materials under discussion are amorphous glass-like 
disordered bodies" for which "Nernst's theorem is not valid^ This remark is erro- 
neous in two aspects. First, as discussed above, there are dielectric materials with 
an ordered structure for which n does not vanish with T (dielectric like semimet- 
als, for instance). Second, the violation of Nernst's theorem for disordered bodies is 
irrelevant to the Casimir entropy discussed here. As is correctly stated in Ref. [THl 
nonzero entropy of glass plates at T = is connected with the fact that they are 
simply not at an equilibrium state at low T. This entropy does not depend on 
separation between the plates and, thus, is not in contradiction with the Nernst 
theorem. By contrast, the Casimir entropy (|3ip is nonzero and depends on separa- 
tion distance (i.e., on the volume of the system) for both ordered and disordered 
materials of plates (provided n does not vanish with T) which is in contradiction 
with the Nernst theorem. 

Thus, the substitution of the modified reflection coefficients into the Lifshitz 
formula leads to contradictions with thermodynamics. The physical reason for this 
is that the Lifshitz theory describes a system in thermal equilibrium whereas the 
modified reflection coefficients were obtained for a system that includes both drift 
and diffusion currents which are irreversible processes out of thermal equilibrium. 
As was emphasized at the beginning of this section, the Boltzmann transport equa- 
tion describes irreversible processes which occur, for instance, in an external electric 
fleld. It is a far reaching extrapolation to apply this equation to the fluctuating elec- 
tromagnetic field, as is done in Ref. [151 In an external electric field the system goes 
out of thermal equilibrium and the fluctuation-dissipation theorem is violated. This 
cannot happen and docs not happen in the presence of fluctuating electromagnetic 
fields. One can conclude that the thermodynamic puzzles discussed above are ar- 
tifacts of the application of the Lifshitz theory and Boltzmann transport equation 
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Fig. 3. The differences of the theoretical and the mean experimental Casimir pressures between 
the two Au plates versus separation are shown as dots. The theoretical results are calculated using 
the modified Lifshitz theory with the inclusion of screening effects. The solid and dashed lines 
indicate the boundaries of 95% and 99.9% confidence intervals, respectively. 




outside of their range of applicability. 



5. Experimental Tests for the Influence of Screening Effects on the 
Casimir Force 

The theoretical predictions made by the approaches including screening effects in 
the Lifshitz theory can be compared with the recent experimental results from the 
measurement of the Casimir force between metal-met al a nd metal-semiconductor 
test bodies. In fact all three theoretical approacheJl^^tl^ discussed above lead to 
almost coincident predictions. This is because the contributions from the zero fre- 
quency term in all these approaches are exactly the same and the contributions 
from all nonzero Matsubara frequencies are approximately equal. 

We begin with the most precise experiment on the indirect measurement of the 
Casimir press ure between two gold plates by means of a micromechanical torsional 
oscillator .1^^ In this experiment, the configuration of a sphere above a plate was 
used in the dynamic regime to determine the equivalent Casimir pressure in the 
configuration of two plates using the proximity force approximation. (This experi- 
ment was already mentioned in Sec. 2 in the comparison between the Drude model 
and the generalized plasma-like model in combination with the standard Lifshitz 
theory.) Note that Fig. 1(a) can be also used to compare the experimental data 
shown as crosses with the predictions of the theoretical approaches taking screening 
effects into account .'^^'"'1^ The point is that at separations larger than 300 nm the 
approaches including the screening effects with k = ktf lead to the same compu- 
tational results for the Casimir pressure as the Drude model approach shown as 
the dark-gray band (the latter is plotted in the region from 500 to 600 nm). As is 
seen in Fig. 1(a), the approaches taking into account screening effects are experi- 
mentally excluded at a 95% confidence level. In Fig. 3 we use another method of 
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comparison of these theoretical approaches with data over the entire measurement 
range from 162 to 750 nm. Here, the differences between the theoretical and the 
mean experimental Casimir pressures are shown as dots. The solid and dashed lines 
indicate the boundaries of 95% and 99.9% confidence intervals, respectively. As is 
seen in the figure, the dots are outside the 95% confidence interval over the entire 
measurement range. At separations from 162 to 640 nm the theoretical approaches 
including the screening effects are experimentally excluded at the 99.9% confidence 
level. Recently it was claimecPSI that the electrostatic calibrations in the experi- 
ments of Refs. [Tj (8, do not take into account "important systematic effects" . This 
claim remains unjustified until some specific objection to the calibration procedures 
of Refs. [71 [8] described in more detail in Refs. |46l [47l is presented. It was claimed 
also^Slthat the comparison of experiment with theory, like in Figs. 1, 3, is irrelevant 
because the optical properties of the Au films used were not measured but taken 
from tables. However, computations performed in Ref. [49] demonstrated that the 
use of any alternative set of optical data only increases disagreement between the 
Drude model approach or approaches of Refs. [TlHIBl and the experimental data. 

Next we compare the theoretical predictions with the inclusion of the screening 
effects with the measurement data of the experiment on the modulation of the 
Casimir force with light.---— (This experiment was also mentioned in Sec. 2.) 
Here, the difference of the Casimir forces F^is between an Au coated sphere and 
a Si plate was measured in the presence, and in the absence, of laser light on the 
plate. In the absence of light, the concentration of charge carriers in Si was much 
below the critical value, i.e.. Si was in a dielectric state. In the presence of light, the 
concentration of charge carriers in Si was above the critical value, i.e.. Si was in a 
metallic state (all values of respective parameters are listed in Ref. 1^ . 

The theoretical predictions with the inclusion of screening effects are numerically 
almost the same, as in the standard Lifshitz theory, for the case when the Si plate 
is in the bright phase, but differ measurably when the Si plate is in the dark phase 
in comparison with the calculation where the conductivity of Si in the dark phase 
is simply neglected (see Sec. 2). It is not possible, however, to conclusively compare 
experiment with theory taking the screening effects into account at a 95% confidence 
level. Below we compare the experimental data whose total experimental errors are 
determined at a 70% confidence level with the theoretical bands computed with 
inclusion of screening effects. Note that in the dark phase k ~ kdh was used in the 
computations. In the bright phase the computational results for the Casimir force 
are almost the same for both k, = kdh and k = ktf- The width of the bands is found 
at the same 70% confidence level as the experimental errors. This width is mostly 
determined by the error in the charge carrier densities used in the computations. 
In Fig. 4(a,b) the experimental difference Casimir forces, F^^^^ (the force in the 
bright phase minus the force in the dark phase) are shown as crosses for different 
measurements with absorbed power = 9.3 mW and 8.5 mW, respectively. The 
theoretical bands for i^^j^"' computed with the inclusion of the screening effects lie 
in between the dashed lines. As is seen in Fig. 4(a,b), the theoretical approach taking 
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Fig. 4. (a) The measured force differences arc shown as crosses versus separation a for the ab- 
sorbed power (a) = 9.3 mW and (b) Pw = 8.5 mW. The theoretical results calculated with 
inclusion of screening effects lie between the dashed lines. 



into account the screening effects is excluded by the data at almost all separation 
distances within the measurement range from 100 to 300 nm. 

In Fig. 5 the experimental data from one more repetition of the same experiment 
with a lower absorbed power {Pw = 4.7 mW) are compared with the theoretical 
approaches taking the screening effects into account. In Fig. 5(a) the experimental 
data are shown as crosses. The theoretical results computed with the inclusion of 
the screening effects belong to the band in between the two dashed lines. As is seen 
in Fig. 5(a), below 200 nm and above 275 nm the theory is inconsistent with data. 
In Fig. 5(b) the same data are compared with the same theoretical approach using 
another method of comparison. Here, theoretical minus experimental differences of 
the Casimir forces are plotted as dots. The solid line shows the upper boundary of 
the 70% confidence intervals taking into account all experimental and theoretical 
errors. It is seen that below 185 nm and above 280 nm all dots are outside the 
confidence intervals. As a result. Figs. 4(a,b) and 5(a,b) conclusively indicate that 
the theoretical approaches taking into account screening effects are excluded by the 
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Fig. 5. (a) The measured force differences are shown as crosses versus separation a for the ab- 
sorbed 4.7 mW. The theoretical results calculated with inclusion of screening effects lie 
between the dashed lines, (b) Theoretical minus mean experimental differences of the Casimir force 
for the Lifshitz theory with inclusion of screening effects. The solid line indicates 70% confidence 
intervals. 
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data of the optical modulation experimeiilP32Hl at a 70% confidence level. 

In contrast to this conclusion, it was claimecfl^ that "the experiment can hardly 
distinguish between nonlocal theory and local theory with zero conductivity of Si." 
(We recall that according to Sec. 2, Fig. 2, the local theory with dc conductivity of Si 
neglected in the dark phase is in perfect agreement with the data.) The comparison 
between the theoretical and experimental results in Rcf. 16 is, however, irregular.^^SI 
In Fig. 1(a) of Ref. [TOI the experimental data [taken from Fig. 1(a) of Ref. [251 with 
no indication of the source] are shown with errors determined at a 70% confidence 
level. In the same figure the theoretical band for the nonlocal approach including 
the screening effects is obtained from the uncertainty in n, An = 0.4 x lO^^cm"'^, 
determinecP^ at a 95% confidence level. The reader of Rcf. 16 is not informed 
about the confidence levels used. This comparison of experiment with theory is 
thus confusing. In our Fig. 4(a) the same data are compared with the predictions of 
the nonlocal approach (the band between the dashed lines), where the band width 
is determined at the same 70% confidence level, i.e., with An = 0.3 x 10^® cm~'^, as 
the errors of the data. It can be clearly seen that the nonlocal approach is excluded 
by the data at a 70% confidence level. 

Regarding Fig. 1(b) in Ref. [121 it shows both the data and the theoretical bands 
at a 95% confidence level and restricts the region of separations only from 100 to 
150 nm. However ,^221 these data cannot be conclusively compared with the nonlocal 
approach at such a high confidence. The correct comparison between the data and 
the nonlocal approach for this measurement set of the optical modulation experi- 
ment over the entire separation range is presented in our Fig. 4(b). It can be clearly 
seen that data are inconsistent with the nonlocal approach including the screening 
effects over a wide range of separations. 

6. Conclusions and Discussion 

In the above, we have reviewed problems with the Lifshitz theory when the drift 
current of conduction electrons in metals, or the dc conductivity at T 7^ in di- 
electrics, are included in the model of the dielectric response. It was shown that 
these phenomena violate thermal equilibrium which is the basic applicability con- 
dition of the Lifshitz theory. Because of this, the inclusion of the drift current into 
the Lifshitz formula results in thermodynamic and experimental inconsistencies. 

The attempts to avoid these problems introduce modified refiection coefficients 
taking into account screening effects and diffusion currents.^l^l"^^ We show (see 
also Refs. l3&lf4T|) that the Lifshitz formula combined with the modified reflection 
coefficients leads to nonzero Casimir entropy at T = depending on the separation 
distance between the plates. Thus, the proposed approaches violate the third law of 
thermodynamics. These approaches are also demonstrated to be inconsistent with 
the measurement data of two experiments. For metal-metal test bodies they are 
excluded experimentally at a 99.9% confidence level. In the configuration of metal- 
semiconductor test bodies the exclusion is confirmed at a 70% confidence level. 
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The reason for the failure of the theoretical approaches, including the screening 
effects and diffusion currents into the Lifshitz theory, is that these processes violate 
thermal equilibrium which is the basic applicability condition of the Lifshitz theory. 
Both the drift and diffusion currents are irreversible phenomena which take place out 
of thermal equilibrium. By contrast, the dispersion forces are physical phenomena of 
fluctuating nature which go on in thermal equilibrium. One arrives at the conclusion 
that there is a deep difference between cixternal and fluctuating electromagnetic 
fields which might not be sufficiently reflected in the mathematical formalism of 
quantum theory. Future studies will show how important this conclusion is for 
wider ranges of physical phenomena beyond the scope of dispersion forces. 
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